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Entanglement properties of two uncoupled atoms embedded in a coherent field distribution 
through one quantum transition process is studied. A case of non-linear Hamiltonian of the problem 
is considered through which the effect of a non-linear media is illustrated. Moreover, the effect of 
the frequency difference between the interatomic transition and the electromagnetic field is also an- 
alyzed. We show that, adjusting the considered parametres of the non-linear media and frequency 
difference leads to a strong control of the degree of entanglement where excellent periodicity of 
entanglement evolution can be obtained which is very important in predicting the behavior of trans- 
mitted information through the application of various information processing schemes. We present a 
' detailed and comparative study of atom-atom entanglement for two cases corresponding to different 

injections of the two atoms into the cavity field. Moreover, we present an answer to the question: 
How does the quantum phase space structure for a composite system relate to the entanglement 
characteristics of the corresponding quantum system? We demonstrate how the entanglement in 
nonlinear tripartite systems can be associated with a derealization in the phase space distribution. 
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I. INTRODUCTION 



The recent rapid development of quantum information theory has largely stimulated research on nonclassical phe- 
nomena, with the main focus on the generation of entangled states that are required for tasks such as quantum 
£f) ' teleportation, dense coding, or certain types of quantum key distribution protocols. In recent years, there are much 
interest to study entanglement of two-level systems. These systems gain such importance because of the efficiency in 
representing information in most quantum information processing schemes. Particular and great interest is devoted 
to the generation of entangled states in two-atom systems, since they can represent two qubits that considered the 
C*~) 1 base blocks of building the quantum gates that are essential to implement various quantum protocols. Furthermore, 
there have been many proposals for generating atomic entanglement and entanglement between cavity modes through 
atom-photon interaction [ll— Tl6j) - Some notable experimental demonstrations have also been performed fl7] . 

However, the physical nature of the interacting objects and the character of their mutual coupling control strongly 
. the degree of quantum entanglement. A number of studies have shown that entanglement can be created between to 
objects which do not interact directly with each other but interact with either a common field or heat bath or thermal 
rS ; cavity field [HI GHH. 

The formation of atom-photon entanglement and the subsequent generation of correlations between spatially sepa- 
rated atoms have been shown using the micromaser [l?], l2ll - |24j . The effectiveness of micromaser [25rl27j in generating 
entanglement has been appreciated as it is strongly considered as a practical device for processing information. It 



stores radiation for times significantly longer than the duration of the interaction with any single atom 27] . The 
interaction of an atom with the intracavity field of a micromaser will, generally, leave the atom-field system in an en- 
tangled state. The long cavity lifetime means that the memory of this entanglement can influence the interaction with 
subsequent atoms and a nonlocal correlations between these successive atoms can be induced leading to a violation 
of Bell's inequality [27] . 

To our knowledge, the study of such a problem when the cavity field contains a coherent distribution is still 
insufficient. The major difficulty may be due to the complex calculations and long analytical mathematics, specially 
when the system Hamiltonian becomes nonlinear. We are interested in the opportunity to use the micromaser with 
coherent distribution as a source of nonlocal entanglement and answer the important questions: how and to what 
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extent can the degree of atom-atom entanglement be cotrolled in these circumstances? Also, can the Husimi Q- 
distribution be used as a helpful method to quantify nonlocal entanglement dynamics? When applied to a system 
of two qubits, its shown that the concurrence H may be expressed in terms of the second moment of the Husimi 
Q-function. More strikingly, when applied to a system of three qubits, the expression for the moment contains all 
three bipartite concurrence terms. In this way this measure captures all classes of entanglement, not only bipartite 
or otherwise [Hj]. 

Our purpose here is to trace analytically and numerically the problem of a nonlocal entanglement created between 
two successive atoms (two-qubits) thr oug h surroun ding the coherent cavity field with a nonlinear medium, namely, 
Kerr-like medium. The Kerr medium [l5l EE [20t l29l - [3lj can be modelled as an anharmonic oscillator with frequency 
uj. Physically this model may be realized as if the cavity contains two different species of atoms, one of which behaves 
like a two- level atom and the other behaves like an anharmonic oscillator in the single-mode field of frequency uj 32] . 
Such a model is interesting by itself. The cavity mode is coupled to the Kerr medium as well as to the two-level 
atoms. A Kerr-like medium can be useful in many respects, such as detection of nonclassical states (33[, quantum 
nondemolition measurement [jjij], investigation of quantum fluctuations [35j], generation of entangled macroscopic 
quantum states [36L l37j. and quantum information processing (38l. l39l|. 

An outline of paper is as follows: In Sec. II we describe in details the physical mechanism of the interaction and 
introduce the modified model of the total system that accommodates processes such as nonlocal interaction in nonlinear 
media. We obtain the wave function of the first part of interaction on which the final wave function of the total system 
depends. In Sec. Ill, we describe the technique we are going to use in quantifying the degree of entanglement. In 
addition, we calculate the the eigenvalues of the non-Hermitian operator that used to compute the mathematical form 
of entanglement measure. Different kinds of atoms injection into the cavity are introduced in sections IV and V, where 
the corresponding wave functions and their associated density matrix operators are computed. In Set. VI, we give 
a little historical review about the Husimi Q-function, where that function is calculated for different circumstances 
separately with detailed discusion included. Our concluding remarks are shown in Set. VII. 



II. THE FULL SYSTEM AND ITS SOLUTION 



We consider two two-level atoms with the same transition frequency traverse a high Q (Q ss 10 9 ) single mode cavity 
one after the other. The flight time of the atoms is, however, long enough so that there is no appreciable overlap 
between the atomic wave functions [HI, 

EE E3, EES]. With this assumption we ignore the possibility of energy 
exchange between the atoms, although, secondary correlations develop between them. The entanglement of their 
wave functions with the cavity photons can be used to formulate local-realist bounds on the detection probabilities 
for the two atoms [13, EH ■ The generation of nonlocal correlations between the two atomic states emerging from the 
cavity can in general be understood using the Horodecki theorem j4(|. It is assumed that the atom field interaction 
time is shorter than the lifetime of the cavity, so that the cavity relaxation will not be considered. The cavity field 
is assumed to be filled with a nonlinear medium, namely, Kerr medium [IE EE HE [29 31]. Also, it is assumed that 
the eigenfrequency of the atomic subsystem differs from that of cavity field subsystem. Under the rotating wave 
approximation (TWA), the Hamiltonian in the interaction picture (assuming that H = 1) 

H mt = y a z + f(x, n) + X (oL a f +d + a), (1) 

where /(x,/i) = x(n 2 — n) represents the nonlinear term with, 

f(x,n) | n) = [x(n 2 - rij] | n) = f(x,n) | n), (2) 

where a (cf) is bosonic annihilation (creation) operator for the single mode field of frequency w and A is atoms-field 
coupling constant, the operators t? + = |+)(— | , <?_ = |— )(+| and a z = |+)(+| — |— )(— | represent, respectively, 
the raising, lowering and population atomic operators which obey the commutation relations [f7 + ,er_] = a z , while 
A = uo — uj is the detuning parameter which represents the difference between the atomic and cavity subsystems 
eig enfrequencies. We denote by x the dispersive part of the third order susceptibility of the Kerr-like medium [l5l EE 

EE El ml 

Since the state vector, \iPf(0)}, of the field is represented by a linear superposition of the number state \n) 

oo 

|<M0)) = £c», (3) 

and the first atom that traverse the cavity is assumed in its upper state | +), i. e., 

IMO)) =1 +> (4) 
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the initial state vector of the interacting first-atom-field system is given by 

oo 

\M0)) = |Vm(0)) ® \M0)) = £ C>, +}, (5) 

where |n) is an eigenstate of the number operator a^a = n; a'a\n) = n\n), and C n is, in general, complex where the 
square of its modulus gives the probability of the coherent field, with mean photon \a\ 2 = n, to have n photons by 
the relation 

f.n 

P(n)H(n|MO))| 2 Hn 2 = e^-. (6) 

n! 

Under the assumption that, first atom traverse the cavity in its upper state | +), the joint state vector of the field and 
the first atom at any instant of time t can be obtained from the solution of the time-dependent Schrodinger equation 

ij t \^A F (t)) =H |VafW). (7) 
The time-dependent wave function of the first-atom-field system takes the form 

\^ AF {t)) =| Ul{t))+ | U n _ +1 {t)) (8) 
which can be simply obtained by recalling the initial condition ([5]) and solving the Schrodinger equation ([7J), where 

oo 

\UZ(t)) = J2C n T 1 (n,t)\n,+), (9) 



I U"l +1 (t)) =J2 C n r 2 (n + l,t) |n + l,-), (10) 

n 

with the amplitudes Ti(n,t) and T 2 (n + l,t) are, respectively, given by 

2 

T 1 {n,t)=Y J F ] {nY H{n) \ (11) 

3=1 

and 

T 2 (n + l,t) = - Y Fj ^W + ^-tofaU (12) 
?=i 



with Fj(n) is given by 



where the angles ^\^{ n ) are given by 



W (n) + «(n) Mi = 1 , 2) (13 ) 



a(?i) + 7(71) / / a(n) — 7(71) x " 



Mi,a(n) = - W 2 " ' ±\l\ " \ — ) +\ 2 (n + l). (14) 



with a(n) and 7(71) read, respectively 



a ( n ) = y + Xn{n - 1), (15) 



7(") = —g + + !)■ (16) 
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As already indicated above, we consider a pair of two-level atoms going through the cavity mode one after another. 
Moreover, under the assumption that time of flight through the cavity t is the same for every atom [H d |H lioMiEj , 
where the joint state vector of both atoms and the field may be denoted by \ifjAAF(t)}, the corresponding atom-atom- 
ficld pure-state density operator may be written as 

PAAF{t) = \TpAAF(t))(lpAAF(t)\, (17) 

where, the joint time-evolved wave vector, \ipAA F (t)) , of the tripartite system, i.e., the two atoms and the cavity field 
after the second atom leaves the cavity is obtained by solving the Schrodinger equation 

lj t \lpAA F {t)) =H mt \i>AA F (t)). (18) 

It is worth to note that within the delay time between the two atoms the field evolves towards a thermal steady state, 
moreover, repetition of the instant in which the later atoms enter the cavity means the same field repeats at this 
instants precisely when successive atoms exit the cavity (25| . 

In order to quantify the degree of entanglement between the two atoms, the field variables must be traced out. One 
may write the reduced mixed-state density matrix of the two atoms after taking the trace over the field variables as 

PAA{t) = Tip [pAA F (t)]. (19) 

Matrix representation of the reduced density operator, Eq. (|19l) becomes 



pAA(t) 



Pll Pl2 Pl3 Pl4 \ 

P21 P22 P23 P24 

P31 P32 P33 P34 

P41 P42 P43 P44 



(20) 



III. ENTANGLEMENT MEASURE 

For bipartite pure states, the partial (von Neumann) entropy of the reduced density matrices can provide a good 
measure of entanglement. However, for mixed states von Neumann entropy fails, because it can not distinguish 
classical and quantum mechanical correlations. For mixed states, the entanglement can be measured as the average 
entanglement of its pure-state decompositions Ep(p) 

E F {p)=min^PiE[iPi), (21) 

i 

with E(ipi) being the entanglement measure for the pure state ipi corresponding to all the possible decompositions 
p = '^2 i Pi\ipi)('4 , i\- The existence of an infinite number of decompositions makes their minimization over this set 
difficult. Wootters Q succeeded in deriving an analytical solution to this difficult minimization procedure in terms 
of the eigenvalues of the non-Hermitian operator 

T = pp, (22) 
where the tilde denotes the spin-flip of the quantum state, which is defined as 

p = {(T x ®C7 x )p*{cr x ®(j x ), (23) 
where a x is the Pauli matrix, and p* is the complex conjugate of p where both are expressed in a fixed basis such as 
In terms of the eigenvalues of T = pp, Ep(p) (known as the entanglement of formation) takes the form 



1 + 



(24) 



E F {p) = h 

where C(p) is called the concurrence and is defined as 

C(p) = max (0, y/8x - y '£* - y/£ 3 - V^4 ) , (25) 
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with the S'a representing the eigenvalues of T = pp in descending order, and, 

Kv) = -y lo s y- ( l - y) 1 °g( 1 - y)> 



(26) 



is the binary entropy. The concurrence is associated with the entanglement of formation Ep(p), Eq. (|24[) . but it is by 
itself a good measure for entanglement. The range of concurrence is from to 1. For unentangled atoms C(p) = 
whereas C{p) = 1 for maximally entangled atoms. 

Recalling f[23]) and assuming that that the matrix T = paa{<Jx ® &x)paa(' j x ® a x), needed for calculation of the 
concurrence, has the form 



/ Tn T12 T13 Ti4 

T21 T 2 2 T23 T 2 4 

T31 T 32 T33 T34 

V T41 T42 T43 T44 



(27) 



the next step is to find the eigenvalues of the above matrix T. To achieve our goal, we are supposed to solve the 
characteristic equation of 



Det(T - £1) = 0, 
Equation (j2"5)) is a polynomial equation of degree 4 

£ i + c 3 £ 3 + c 2 £ 2 + ci£ + c = 0, 

with 

C3 = —Tii — T 2 2 — T33 — ^44, 



(28) 
(29) 
(30) 



c 2 = -|T\ 3 | 2 - |Ti4| 2 - |Ti 2 | 2 - |T 34 | 2 - |T 24 | 2 + (Tn + T 22 )(T 33 + T 44 ) + TnT 22 + T33T44, 



(31) 



ci = |Ti 2 | 2 (T 33 + T 44 ) + |T 23 | 2 (Tn + T 4 4)|T 24 | 2 (Tii + T 33 ) + |T 34 | 2 (Tii + T 22 )|Ti 3 | 2 (T 22 + T 44 ) + |Ti 4 | 2 (T 22 + T 33 ) 
— 3^(T 2 3T 3 4T4 2 ) — ^(T 2 iT32Ti3) — 5R(T 2 iT4 2 Ti4) — 3fJ(T3iTi4T4 3 ) — Tu(TnT3 3 + T 22 T 33 + TnT 22 ), (32) 
co = TnT 2 2T 33 T 44 + |Ti 3 | 2 |T 24 | 2 + |Ti 4 | 2 |T 23 | 2 - TiiT 22 |T 34 | 2 - TiiT 4 4|T 23 | 2 - T 33 T 4i \T 12 \ 2 
— TiiT 33 |T 2 4| 2 — T 22 T33|Ti4| 2 + Tn3?(T 32 T 2 4T43) + T 22 3?(T 3 iTi4T4 3 ) + T 33 3fi(T 2 iT4 2 Ti4) + T i4 ^.{T 3 2Ti 3 T 43 ) 



The roots of Eq. 



— 5R(T 2 iT 32 Ti4T 4 3) — 3?(T 2 iT42Ti 3 T34) — Sft^iT^TuT^). 
are as follows 



£1 



C3 



40 



(33) 



(34) 



C3 



(35) 



Sa = - 



c;s 



AN 



(36) 



* = " 4 



C3 



W 

T + 4V 



(37) 



with 



where 



and 



with 



and 



where 



o.i/ft+ft+a (38) 



N ^q 1± o i -0 1 m 



Ox = 3c§ - 8c 2 - V 1/3 (40) 
2 = 4V- 1 / 3 [3ciC3 - 12c - c|] (41) 



03 = y^[4c 3 c 2 -ci-8 Cl ] (42) 



Vi = 8c| + 36[3( C 2 + cocj) - Cl c 2 c 3 - c c 2 ] (44) 

V 2 = -Mc 2 [c QCl cl + cjc 3 + 8c (c? + c cl)} (45) 

Vs = 3 Cl C3[2co(4 C 2 + 3cic 3 ) - ci^Cs] (46) 

V 4 = 12[( Cl c 3 ) 3 + 4c oC 2(8c - c§) + c£(c? + cqc§)] (47) 

V 5 = cl[81cl + 576c lC3 - 768c c ] (48) 

W = W 1 +W 2 (49) 

Wi = 3c 2 3 - 8c 2 + 2V 1/3 (50) 

W 2 =8V- 1 ^[c 2 2 + 3(4c -c 1 c 3 )} (51) 

In the following we consider two different cases of injection of the second atom into the cavity subsystem that give 
us clear insight into the mechanism of how and to what extent the degree of entanglement can be controlled and 
enhanced. 
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IV. INJECTION OF TWO EXCITED ATOMS SUCCESSFULLY 



In this case, the time dependent state vector of the full system can be expressed in the form 



\ipAAF(t)) = 



C/ + (t))|n> + 



U ± (t))\n + l) 



U ^ (f))|n + l> + 



U_ (t)) \n + 2), 



which can be simply obtained by solving the Shrodingcr equation (jl8)) , where 

u^(t)\=J2c n [iMm)] 2 | + ), 



(52) 



(53) 



U ± (*)> =E C " ^(Mj^fr + M)! ± >, 



(54) 



C/ T (*)) =E C " ri(n+l,t)r 2 (n + l,t)| T ), 



(55) 



C/ = (t)) =E C " r 2 (n+l,t)r 2 (n + 2,t)| = ), 



(56) 



Using (fT7| to obtain the full density operator and after the application of (Q~9|) to obtain the reduced atomic density 
operator which can be put, after using the notations | =|= ) =| 1), | ± ) =| 2), | =p ) =| 3), | = ) =| 4), in the form 



where 



pAA(t) = #j (*)!*> Oil 



p u (t) = ^|C7"| 2 |r 1 (n,t)| 4 , 



(57) 



(58) 



M*) = ^C n+1 C*"[r 1 (n + l,i)] 2 r i (n,t)r*(n + 1,<), 



(59) 



Pis(t) = E C n+1 C* n |r!(n + l,i)| 2 ri(n + l,t)T* 2 (n + l,t), 



(60) 



P14 (*) = E C" +2 C*" [r x (n + 2, t)] 2 r*(n + 1, t)r*(n + 2, t), 



(61) 



P2 2 (t) = ^|c"| 2 |r 1 (n,t)| 2 |r 2 (n + i,t)| 



(62) 



P23 (t) = 1 C 1 2 1 r 2 (" + 1 . l ) T r i OH (n + 1 , t) , 



(63) 



p 24 (<) = E ^'^T^n + 1, t)r5(n + 1, i)|r 2 (n + 2, t)| 



(64) 



pasCf) = El c '1 2 | r i(" + M)| 2 |r 2 (n + i,t)\ 2 , 
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(65) 



P3*(*) = 22 C n+2 C* n Ti{n + 1, t)T* 2 {n + 1, t)|r 2 (n + 2, t)| , (66) 



P44 



(t) = 1 C " 1 2 1 r 2 (« + 1 . *) 1 2 1 r 2 (« + 2, t) 1 2 , (67) 



For x = A = 0.0, the results reported in Ref. are straightforward obtained, which analyzed in details but for low 
and high average photon number fi. 

In order to discuss the dynamical feasibility of our scheme in the micromaser regime, in which the lifetime of the 




FIG. 1: The concurrence C(\t/n) (solid red curve) and the total population p22(Ai/7r) + p33(\t/n) (dashed blue curve) for 
ri = 10, S = 0.0 where (a) xA = 0.0 (b) X /A = 0.2 (c) xA = 0.5 (c) xA = 1-0 

cavity is longer than atom held interaction time, the concurrence C(\t/ir) and the corresponding total population 
/922(A£/7r) + ,033 ( At A) > as a reference state, are depicted simultaneously in figures 1, 2 and 3 for the Rabi angle 
Xt/ir £ [0, 16]. The mean photon number h is taken about n — 10, the detuning parameter S = A/A carries various 
values such as S = 0.0,0.5,1,10,15 and 20, while the Kerr parameter x/A varies such as x/A = 0.0,0.2,0.5 and 
x/A = 1.0. From Fig. la, where the interaction occurs resonantly in the absence of the nonlinear media, one can 
see clearly that the concurrence C exhibits irregular collections of oscillations where the total population oscillates 
but with amplitudes maximize (no longer exceed the value 0.25) when P22 + P33 ~ 0.4; the atoms in the mixed state 
I ± ) + I =p )■ Thus, the lack of the entanglement between the two atoms can be attributed to a large population of 
the product states | ± ) and | =p ) which, in turn, decay on the same time scale as C . For significantly high average 
photon number n, it was shown that [111 ] . quantum effects which are predominant primarily when the photon number 
is low, help to increase the peak value of the concurrence C. 
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FIG. 2: The same as Fig. 1 but for xA = 0.0 where (a) S = 0.5 (b) 5 = 1.0 (c) S = 10.0 (d) 5 = 15.0 (e) S = 20 



On resonance interaction between the atoms and a cavity filled with a Kerr medium of a parameter of order 10 _1 A, 
such as x — 0.2A, the total population shows pseudo periodical behavior with amplitudes minima decrease as time 
evolves where the revival period, ta, elongates as time goes on. The concurrence C loses its irregular collections where 
instant oscillations appear after wide time intervals accompany with amplitudes decrease (sa 0.3 in the beginning) as 
evolves , Fig. lb. A significant and interesting distinction in the case when \ = 0-5A is that both C and P22 + P33 
display a perfect ir periodicity of oscillations. The total population P22 + P33 exhibits collapses and revivals with the 
revival period is shorter and each revival is a perfect carbon copy of its predecessor. The concurrence C shows 
oscillation packets composed of three peaks with the strongest, = 2.5, lies at the center where the total population 
collapses. Moreover, the concurrence C minimizes where the total populations P22 + P33 revives while the strong 
entanglement and collapse of population occur in parallel manner, see Fig. lc. This periodicity is due to the fact that 
analytical expressions of both C and P22+/O33 can be summed exactly. Thus, the matching value of the Kerr parameter 
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FIG. 3: The same as Fig. 2 but when x = 0.5 and (a) 5 = 1.0 (b) 6 = 10.0 



plays more efficient rule than that of quantum effects that leads to increase the peak value of the concurrence C [TTJj . 
Opposite results are obtained when a considerable strength of the Kerr medium is applied to the cavity field. In this 
case, a pseudo periodicity for C is observed where the total population preserves its periodicity accompanied with 
clear decreasing in its amplitudes. In addition, one can't build a clear insight about the relation between curves of 
both C and P22 + P33, Fig. Id. Thus, from the above it is found that the Kerr medium can be used to adjust the 
entanglement degree between qubits, which may provide some reference and even basis for the quantum control of 
multiple qubits. 

A surprise is that the case of a non-resonance interaction does not destroy the entanglement phenomenon, see Fig. 
2. For small detuning such that 6 = 0.5, there is no significant change, Fig. 2a, but when S increases to unity, the 
correlations time between the atoms becomes longer with the gradual appearance of nonzero values of C, see Fig. 2b. 
Strictly speaking, the destruction of entanglement occurs for the case when the detuning parameter S is absent but 
not the Kerr parameter x/A, Fig. Id. 

The most intersting and surprising is the case when the inter-atomic transition is considered to be significantly 
high compared with the cavity frequency such that 5 — 10.0, 15.0 and 20.0 where the nonlinear media is considered 
to be absent. This is revealed in Figs. 2c-e. We see that the total population exhibits the known collape and revival 
phenomena periodically. Also, we notice that for the increase of detuning parameter 6, the collapses become longer 
and the beginning strong overlap of the revivals becomes weaker; while the amplitude of the revivals decreases as 
values of 6 increase. The entanglement evolution divided into three regions. The first one is at the half of the 
collapse time, where the entanglement degree oscillates with bigger amplitude for short time scale. The oscillation 
maximum is related crucially to the value of <5, where it occurs at time Xt = (| + l)7r. The last two regions 
are in the beginning and end of the collape time. In this case, we see that the atom-atom entanglement remain 
close to its maximum and show intrinsic oscillations only. This is in fact is true for all values of the detuning 
parameter. The detuning between the field and the atoms leads to an emergence of collapses and revivals, this 
statement is always true only in the absence of the nonlinear media, see Figs. 2c-e. Overall, we find that the 
detuning parameter acts as a control parameter for the atom-atom entanglement, where a perfect periodicity can 
be obtained for suitable high detuning. The revival of the entanglement at time Xt = (| + l)n is more easily 
understood by reference to the density matrix of the system. As expected, there is no entanglement before the 
populations of the state | ± ^ + | zp ^ and the state | ) + | = ) depopulates, but some entanglement appears 

for longer times, Figs. 2b-e, and the concurrence becomes equal to the population of the state | ± ) + 1 =p ) , sec Fig. lc. 



V. INJECTION OF A GROUND ATOM AFTER THE EXCITED ONE 

In this case, the time dependent state vector of the full system can be expressed in the form 



\ipAAF(t)) = 



S ± (t))\n) + 



5 + (f)) |n-l> + 



S _ (t) ) I n + 1} + 



S T (t))\n), 



(68) 
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where 



5 ± (*)) = E c '" r i( ri '*)n("-i ) *)l ± > ) 



s^(t)) = J2 cn ri(n,t)r 2 (n,t)| + 



S_ (t)\ =J2 C " ri(n + l,t)r 2 (n + l,i)| = ), 

/ 7-, 



S T (t))=J2 Cn [T2(n + l,t)] a |T>, 



Following the same procedure as in the previous section, and after using the notations | ± ^ 
| = ) = | 3) , | =p ) = | 4) , the elements of the reduced atomic density matrix ([57]) read 



(t) = Y,\c n \ 2 \r 1 (n,t)\ 2 \r 1 (n~i,t)f 



(69) 
(70) 
(71) 
(72) 

=1 1), | *) =1 2), 
(73) 



Pia(t) = Yl C"C* n+1 ri(n, t)r;(n - 1, t)rj(n + 1, t)r3(n + 1, t), 



(74) 



p 13 (t) = C n+k C* n \T x {n, t^T^n + 1, i)r*(n + 1, t), 



(75) 



p 14 (i) = V |C7 n | 2 r 1 (n,t)rl(n- l,t)[r*(n + l,t)] 2 , 



(76) 



p 22 (t)=E|C"| 2 |r 1 (n,t)| 2 |r 2 (n,t)| 2 , 



(77) 



ft» 3 (t) = E C m+2 C"T 1 (rc + 2, t)r 2 (n + 2, t)r x (n, t)r*(n + 1, t), 



(78) 



p 24 (f) = ]T c™+ 2 c*"|r 2 (n + 1, t)| 2 r x (n + 1, t)rs(n + 1, t), 



(79) 



P33W-El C "| 2 | r i("'*)| 2 | r 2^ + 1 ' < )| 1 



(80) 



p M (t)=J2C n C* n+1 r* 1 (n,t)T 2 (n + l,t)[T* 2 (n + 2,t)} 2 , 



(81) 



P44 (t) = ^|C"| 2 |r 2 (^ + l,t)| 4 , (82) 
To better understand the situation and gain clearer insight into the entanglement dynamics, hereafter we shall analyze 
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FIG. 4: The concurrence C(\t/ir) (solid red curve) and the total population p22(Ai/-7r) + p33(Ai/7r) (dashed blue curve) for 
n = 10, 6 = 0.0 where (a) x/A = 0.0 (b) x/A = 0.2 (c) x/A = 0.5 (c) x/A = 1.0 

the evolution of the concurrence C(Af/7r) and the total populations p22(Xt/ir)+p33(\t/Tr) where the parametres remain 
as the same in the past section. 

Generally, a similar conclusion is reached for the present case but with considerable distinctions in its details. 
As in the case of two excited atoms, there is no initial entanglement between the atoms, and at early times the 
entanglement builds up rapidly to a maximum appearing at short time Xt « 2.5-7T, Fig. 4a. However, comparing 
with the case of two successive excited atoms, shown in the previous section, we see that the entanglement maxima 
that obtained in the present case are generally lower than that obtained with two successive excited atoms, see 
Figs. 1-3 and 4-6. It is interesting to note that after passing through the maximum, the entanglement decays with 
two different time scales, Figs. 2 and 5. This effect is more pronounced if we choose the interaction with atoms to 
be initially in different quantum states. Thus, atom-atom entanglement depends crucially on the coherent transfer 
of the population between the atoms states, specially in the presence of a nonlinear media as well as the detuning 
parameter 5. In this process, the population is efficiently transferred from the more populated state, | 4^ ) + | = ) to 

the state, ± ) + I =p } before it decays to the initial state leading to the enhancement of the entanglement. 



VI. PHASE SPASE DISTRIBUTION 



We devote this section to concentrate on a function, which considered an important one of the functions, W, (Husimi) 
Q, and (Glauber-Sudershan) P functions, that is the Husimi Q-function which has the nice property of being always 
positive and further advantage of being readily measurable by quantum tomographic techniques 51 , [52 1 . It has been 
shown from earlier studies [47H50lj that these quasi-probability functions are important for the statistical description 
of a microscopic system and provide insight into the nonclassical features of the radiation fields. The homodync 



13 




0.4 



0.3 



0.2 
0.1 








6 8 10 12 14 16 

A,t/jt 



FIG. 5: The same as Fig. 4 but for xA = 0.0 where (a) <5 = 0.5 (b) S = 1.0 (c) 5 = 10.0 (d) S = 15.0 (e) 6 = 20 



measurement of an electromagnetic field gives all possible linear combinations of the field quadratures. The average of 
the random outcomes of the measurement is connected with the marginal distribution of any quasi-probability used 
in quantum optics. In fact, Husimi Q-function is not only a convenient tool to calculate the expectation values of 
anti-normally ordered products of operators, but also to give some insight into the mechanism of interaction for the 
model under consideration. The relation between the phase-space measurement; Husimi Q-function; and the classical 
information-theoretic entropy associated with quantum fields was introduced by Wehrl [53j, which referred to as the 
Shannon information of the Husimi Q-function. Thus, Husimi Q-function can be related to quantum entanglement in 
different approaches [53j-|60(. Furthermore, separable state is represented by a localized wave packet in phase space. 
Since coherent states are the most localized states in the Husimi representation, it is argued that [28} derealization 
of the Husimi distribution implies correlation-hence entanglement- between system particles. 
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FIG. 6: The same as Fig. 5 but when x = 0.5 and (a) S = 1.0 (b) S = 10.0 



The Husimi Q-function can be given in the form as [501 , 16 ll |6 



Q(a) 



{a | p F I a) 



(83) 



where pf is the reduced density operator of the cavity field given by tracing over the atomic variables of the full 
density operator (IT7l) . The state | a) represents the well-known coherent state with amplitude a = X + iY. Inserting 
the obtained pp into Eq. (|83|). we can easily obtain the Husimi Q-function of the cavity field 



-((a, + | p F | a, =^ ) + (a, = \ p F \a, =)), 

7T 



where 



|a>=e-l-l*/ 2 £^= 

n V ™ 

n— 



(84) 



(85) 



A. Injection of two excited atoms one by one 

In this case, the Q-function is given by 



Q = 



l 



a, =|= 



U, ,n 



a, 



U ,n + l 



where for the states | =t= ) an( l | = ), Pf reads 



U ^ (t),n)(U^ (t),n 



C/_ (t),ri + 1 )( Z7_ (t),ri + 1 



in this case 



a, =)= 



U (t),n\ =e -(H 2 +l«ol 2 )/2^ 



^[r l(M) ] 2 , 



(86) 



(87) 



(88) 



r in. n f i ) = e -(i«i 2 +i«oi 2 )/ 2 £ i^r 2 (n + 1, t)r a ( B + 2, t), 



(89) 
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FIG. 7: Conour plots for Husimi Q-function for n = 10, 5 = 0.0, \t = 5tt/3 where (a) x/A = 0.0 (b) x/A = 0.2 (c) x/A = 0.5 
(d) X /A = 1.0 
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FIG. 8: The same as Fig. 7 but when 5 = 10.0 and (a) x/A = 0.0 (b) x/A = 1.0 



a = X + iY. (90) 

Here, we present a detailed analytical discussion about how the Husimi Q-function plays an important role in 
quantifying atom-atom entanglement. For this reason, we have pictured the numerical results obtained as contour 
plots in figures 7-9. The parameters was considered as follows: the average photon number n is kept fixed for all plots 
such as n = 10, the detuning parameter carries values such as 5 = 0.0, 10.0, and the Kerr parameter x/A varies as 
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FIG. 9: The same as Fig. 7 but for x/A = 0.5 where (a) At = tt/6 (b) At = tt/4 (c) At = tt/3 (d ) At = vr/2 

x/A = 0.0, 0.2, 0.5, 1.0, where the Rabi angle At = 57r/3 is kept fixed for figures 7 and 8, while for figure 9 it changes 
as At = 7r/6, 7t/4, tt/3, tt/2, where the Kerr parameter x/A is half integer. 

Form figure 7, by which we examine the effect of the Kerr parameter x/A at a fixed Rabi angle At = 57r/3, we 
observe that, in the absence of the Kerr parameter, Fig. 7a, the Husimi Q-function composed of three components, 
two blobs centered at the points (0, ±4) and a symmetric circular peak (Poisson band) centered at the point (4, 0) 
and strongly localized in the right half of the available phase space, i.e., the interaction region [63j], once more, the 
squeezing effect on the number of photons can be observed from the figure. Experimental and theoretical studies 
showed that squeezing of matter wave fields (squeezed states) is closely related to entanglement [64l - [70l |. Moreover, 
we observe that the Husimi distribution blobs delocalized almost around the right half of the phase space, i. e., the 
interaction region of phase space and have structures with almost positive density throughout the available phase 
space. It is worth to note that the bifurcating of the blobs corresponds to the collapse. This can be clearly seen if we 
compare the numerical calculations for total populations and Q-function. It was shown that the bifurcation of the 
Husimi distribution, which is the signature of the formation of Schrodinger cat states, implies correlation, and hence 
entanglement [H, [H, |H, l7ll - [74j |. see also Fig. la. As soon as the effect of the Kerr medium is considered, a clear 
difference in the figures shape appears, where the sensibility of the Husimi Q-function with the nonlinearity parameter 
x/A is very intersting. We observe that the Husimi distribution is taken gradually towards the non-interaction region, 
namely, left half of the available phase space, where we observe that the two blobs rotating in the complex coherent 
state parameter plane in the counterclockwise directions with the same speed where two overlapping regions appear 
with both terms are appreciably vanishing. But note also that the intersection areas with this prescription are not 
in the same position. Upon the increase of the Kerr parameter then one can observe that the overlapping becomes 
dominant and the Husime distribution are taken completely into the interaction region. As a consequence, strong 
entanglement degree between the two atoms can be built up at that time, see Figs, lc and 7c. The conclusion will 
be extremely interesting when we look at the behaviors of both the concurrence C and the Husimi distribution once 
the Kerr medium is taken into account, specially at the moment when the Kerr medium becomes strong at which 
entanglement diminishes and the Husimi distribution is strongly localized and taken into non-interaction region again, 
Figs. Id and 7d. It is important to note that, localization of the Husimi distribution around a fixed point implies low 
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degree of entanglement [H, Iz3| • 

Our observations will become clear when we look at the behaviors of both the concurrence C and the Husimi 
distribution at the moment when detuning parameter possesses a value. The results showed in Figs. 2c, 3b and 8. A 
clear derealization of the Husimi distribution but into the right half of phase space, the interaction region, corresponds 
to longer time interval of entanglement, Figs. 2c and 8a, while the opposite is true, Figs. 3b and 8b. An extremely 
interesting observations can be obtained when we consider various values of the Rabi angle Xt, while the nonlinear 
media is switched on. The results show strongly the role that the Hisimi distribution play in quantifying entanglement, 
where at all moments but not At = n/2, once the suitable Kerr parameter is chosen, stronger entanglement can be 
built up since the Husimi distribution is delocalized but completely spreads between two or three points into the right 
half of phase space, the interaction region, see also Fig. lc. 



B. Injection of excited and ground atoms one by one 



In this case, for the states | 4 1 ) ancl I = ), Pf reads 



PF 



The Q-function is given by 



S^(t),n-l)(S^(t),n-l 



Q 



a, * 



S i , n — 1 



+ 


S = (t), 


n + l\(s_ 


(t),n + l 


2 








+ 


<- 


S _ ,n + lj 


y 











where 



(a, * S ■ (t), n - l) = e -(H 2 +l°Ql 2 )/ 2 ^u^V^ + X ' ') r »(" + *> «) 



(91) 



(92) 



(93) 



a, 



S (*), n + 1 ) = e-d«l 2 +l-l 2 )/ 2 7= ] =n(n, t)T 2 {n + 1, t), (94) 



To obtain clearer insight about the efheiancy of the Husimi Q-function in quantifying entanglement, a comparative 
discusion will be important. Here we examine the behavior of the Husimi function under the same circumstances of 
parameters as in the previous section. This is illustrated in figures 10-12. Generally, the Husimi Q-function developed 
in manner similar to that observed previousley with slight change such that the localized Poisson band barely appear 
which means that nonclassical effect is not pronounced. In this case the entanglement between the two atoms doesn't 
live much more in opposite to the case of two successive excited atoms, where entanglement does appear at almost 
the whole time scale, see Figs, la, 4a, 7a, 10a. The case when we turn on the nonlinear media is also intersting. 
The Husimi distribution moves quickly towards the non-interaction region once the a weak Kerr medium is switched 
on, while a strong localized Poisson band appear clearly in the interaction region which reflects small probability of 
interaction and so lesser times of entanglement, see Figs. 4b, 10b. The suitable choice of the Kerr parameter push 
the Husimi distribution towards the interaction region again, which means stronger entanglement appears at longer 
times, Figs. 4c, 10c. There is no something interesting in case when the nonlinear media is strong, where almost the 
same behavior was seen before, for comparison see Figs. Id, 4d, 7d, lOd. 

However, the effect of the detuning parameter S is quite differente, specially in the presence of strong nonlinear 
media. In such a case, the Husimi distribution displays a strong derealization and spreads between two points in the 
right half of the phase space in case of zero- value nonlinear media accompanied with a significant density throughout 
the available phase space. Moreover, overlapp takes place and the squeezing effect on the number of photons occurs can 
be observed. The strong Kerr media have intersting effect where the Husimi distribution shows ring-like shape around 
the center of the phase space but nevertheless have a significantly different densities throughout the available phase 
space, i.e., interaction and non-interaction regions. The effect of the matching choice the nonlinear media parameter 
x/A on the Husimi distribution at different values Rabi angle Xt points out that enhancement of entanglement depends 
crucially on this suitable choice, see Figs. 4a and 12a-c. 
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FIG. 10: Conour plots for Husimi Q-function for n = 10, 5 = 0.0, At = 5tt/3 where (a) x/A = 0.0 (b) x/A = 0.2 (c) x/A = 0.5 
(d) x/A = 1.0 
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FIG. 11: The same as Fig. 10 but when 8 = 10.0 and (a) x/A = 0.0 (b) x/A = 1.0 

VII. CONCLUSIONS 

We have considered the nature of the entanglement of output two successive atoms from a micromaser cavity for 
pure states input of the atoms when the cavity distribution is coherent. We have considered two different successive 
injections of the initial states of the atoms that traverse the cavity, namely, initially excited atoms and initially excited 
atom followed ground one. We have examined the effect of the frequency difference between the inter-atomic frequency 
and cavity frequency individually and in coexistence of nonlinear media on nonlocal atomic entanglement. We have 
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FIG. 12: The same as Fig. 10 but for x/A = 0.5 where (a) At = tt/6 (b) At = tt/4 (c) At = tt/3 (d ) At = tt/2 



also investigate the role of the Husimi Q-distribution that plays to give clearer insight about entanglement dynamics. 
Our conclusions are Summarized in the following: 

(i) The physical nature of the interacting objects and the character of their mutual coupling control strongly the 
degree of quantum atomic entanglement, in other words, the interaction of a cavity field with two successive excited 
atoms plays more efficient role in producing atomic entanglement than that of interaction with two successive atoms 
traverse the cavity in dissimilar initial states. 

(ii) The nonlinear medium plays an important role in producing atomic entanglement depending on the appropriate 
choice of its parameter. In this case an excellent periodical entanglement can be produced with increased maximum 
value comparing with the case when the nonlinear medium parameter is not appropriate or zero. 

(iii) The extremely interesting is the role of the detuning parameter on producing nonlocal atomic entanglement. 
When the detuning parameter is considered to be considerably high, periodical, long lived entanglement can be 
obtained. 

(iv) Husimi Q-distribution plays a very clear and efficient role in quantifying nonlocal atomic entanglement, where 
bifurcating of Husimi Q-distribution in whole phase space and derealization in the interaction region of the phase 
space corresponding to correlations between particles and as a consequence nolocal atomic entanglement. 



[1] N. A. Peters, T.-C. Wei, and P. G. Kwiat, Phys. Rev. A 70, (2004) 052309. 

[2] C. H. Bennett, H. J. Bernstein, S. Popescu, and B. Schumacher, Phys. Rev. A53, (1996) 2046. 

[3] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. Wootters, Phys. Rev. A 54, (1996) 3824. 

[4] V. Vedral and M. B. Plenio, Phys. Rev. A 57, (1996) 1619. 

[5] V. Vedral, Rev. Mod. Phys. 74, (2002) 197. 

[6] G. Vidal and R. Tarrach, Phys. Rev. A 59, (1999) 141. 

[7] V. Vedral, M. B. Pienio, M. A. Rippin and P. L. Knight, Phys. Rev. Lett. 78, (1995) 2275. 
[8] W. K. Wootters, Phys. Rev. Lett. 80, (1998) 2245. 



20 



[9] B. Ghosh, A. S. Majumdar and N. Nayak, Phys. Rev. A 74, (2006) 052315. 
[10] J. I. Cirac and P. Zoller, Phys. Rev. A 50, (1994) R2799; J. I. Cirac, P. Zoller, H. J. Kimble, and H. Mabuchi, Phys. Rev. 
Lett. 78, (1997) 3221; L.-M. Duan, M. D. Lukin, J. I. Cirac, and P. Zoller, Nature London 414, (2001) 413; E. Solano, G. 
S. Agarwal, and H. Walther, Phys. Rev. Lett. 90, (2003) 027903; S. G. Clark and A. S. Parkins, ibid. 90, (2003) 047905; 
L.-M. Duan, B. Wang, and H. J. Kimble, Phys. Rev. A 72, (2005) 032333; G. S. Agarwal and K. T. Kapale, ibid. 73, 
(2006) 022315. 

[11] B. Ghosh, A. S. Majumdar and N. Nayak, Int. J. Quant. Inf. 5, (2007) 169. 
[12] S. J. D. Phoenix and P. L. Knight, Phys. Rev. A 44, (1991) 6023. 
[13] S. J. D. Phoenix and P. L. Knight, Phys. Rev. Lett. 66, (1991) 2833. 

[14] I. K. Kudryavtsev, A. Lambrecht, H. Moya-Cessa and P. L. knight, J. Mod. Opt, 40, (1993) 1605. 
[15] M. S. Ateto, Int. J. Theor. Phys. 48, (2009) 545 
[16] M. S. Ateto, Int. J. Quant. Inf., 5(4), (2007) 535. 

[17] E. Hagley, X. Maitre, G. Nogues, C. Wunderlich, M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. 79, (1997) 
1; A. Rauschenbeutel, P. Bertet, S. Osnaghi, G. Nogues, M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. A 64, 
(2001) 050301 R; A. Auffeves, P. Maioli, T. Meunier, S. Gleyzes, G. Nogues, M. Brune, J. M. Raimond, and S. Haroche, 
Phys. Rev. Lett. 91, (2003) 230405; S. Nussmann, M. Hijlkema, B. Weber, F. Rohde, G. Rempe, and A. Kuhn, ibid. 95, 

(2005) 173602; D. N. Matsukevich, T. Chaneliere, S. D. Jenkins, S.-Y. Lan, T. A. B. Kennedy, and A. Kuzmich, ibid. 96, 

(2006) 030405. 

[18] M. S. Kim, J. Lee, D. Ahn and P. L. Knight, Phys. Rev. A 65, (2002) 040101(R); L. Zhou, H. S. Song and C. Li, I. Opt. 
B: Quantum Semiclass. Opt. 4, (2002) 425. 
D. Braun, Phys. Rev. Lett. 89, (2002) 277901. 

M. S. Ateto, Applied Mathematics & Information Sciences, 3, (2009) 41. 
M. Loffler, B.-G. Englert and H Walther, Appl. Phys. B 63, (1996) 511. 
A. S. Majumdar and N. Nayak, Phys. Rev. A 64, (2001) 013821. 

G. Rempe, F. Schmidt-Kaler, and H. Walther, Phys. Rev. Lett. 64, (1990) 2783; G. Rempe, and H. Waklther, Phys. Rev. 
A 42, (1990) 1650; H. Paul, and Th. Richter, Optics Commun. 85, (1991) 508; J. D. Cresser, Phys. Rev. A 46, (1992) 
5913. 

P. Masiak, Phys. Rev. A 66, (2002) 023804. 

P. Filipowicz, J. Javanainen, and P. Meystre, Phys. Rev. A 34, (1986) 3077. 
S. J. D. Phoenix and P. L. Knight, Ann. Phys. 186, (1988) 381. 
S. J. D. Phoenix and S. M. Barentt, I. Mod. Opt, 40, (1993) 979. 
A. Sugita, J. Phys. A: Math. Gen. 36, (2003) 9081. 
M.-F. Fang and H.-E. Liu, Phys. Lett. A 200, (1995) 250. 
M. S. Ateto, Int. J. Theor. Phys. 49, (2010) 276. 
C.-Z. Wang, C.-X. Li and G.-C. Guo, Eur. Phys. I. 37, (2006) 267. 
A. Joshi and R. R. Puri, Phys. Rev. A 45, (1992) 5056. 
M. Hillery, Phys. Rev. A 44, (1991) 4578. 

A.N. Chaba, M.J. Collett, D.F. Walls, Phys. Rev. A 46, (1992) 1499. 

R. Zambrini, M. Hoyuelos, A. Gatti, P. Colet, L. Lugiato, M.S. Miguel, Phys. Rev. A 62, (2000) 063801. 

C. C. Gerry, Phys. Rev. A 59, (1999) 4095. 

M. Paternostro, M.S. Kim, B.S. Ham, Phys. Rev. A 67, (2003) 023811. 
J. Pachos, M. Chountasis, Phys. Rev. A 62, (2000) 052318. 

D. Vitali, M. Fortunato, P. Tombesi, Phys. Rev. Lett. 85, (2000) 445. 
A. Datta, B. Ghosh, A. S. Majumdar and N. Nayak, Europhys. Lett. 67, (2004) 934. 
M. Loffler, B.-G. Englert and H Walther, Appl. Phys. B 63, (1996) 511. 

A. S. Majumdar and N. Nayak, Phys. Rev. A 64, (2001) 013821. 

B. -G. Englert T. Gantsog, A. Schenzle, C. Wagner and H Walther, Phys. Rev. A 53, (1996) 4386. 

X. Maitre, E. Hagley, G. Nogues, C. Wunderlich, P. Goy, M. Brune, J. M. Taimond and S. Haroche, Phys. Rev. Lett. 79, 
(1997) 769. 

B. -G. Englert, P. Lougovski, E. Solano and H. Walther, e-print: |quant-ph/02 09128vl 24Sep 2002. 
T. Horodecki, P. Horodecki and M. Horodecki, Phys. Lett. A 200, (1995) 340. 

E. Wigner, Phys. Rev. 40, (1932) 749. 
Z. Wigner, Phys. Chem. B 19, (1932) 203. 
K. E. Cahill and R. J. Glauber, Phys. Rev. 177, (1969) 1882. 

M. Hillary, R. F. O'Conell, M. O. Scully and E. wigner, Phys. Rep. 106, (1984) 121. 
E. L. Bolda, S. M. Tan and D. Walls, Phys. Rev. A 57, (1998) 4686. 
S. Mancini and P. Tombesi, Europhys. Lett. 40, (1997) 352. 
A. Wehrl, Rep. Math. Phys. 16, (1979) 353. 

V. Perinovd, J. Krepelka, J. Perina, J. Luks and P. Szlachetkta, Opt. Acta 33, (1986) 15. 
Hong-yi Fan and Qin Guo, |quant-ph/0611206y i (2006). 

C. P.-Campos, J. R. G.-Alonso, O. Castanos and R. L.-Pena, cond-mat.guant-gas/0910.3256vl (2009). 
Li-yun Hu and Hong-yi Fan, quant-ph/0911.0125vl, (2009). 

A. Miranowicz,J. Bajer, M . R. B. Wahiddin and N. Imoto, J. Phys A: Math. Gen. 34, (2001) 3887. 
A. Miranowicz, H. Matsueda and M. R. B. Wahiddin, I. Phys. A: Math. Gen. 33, (2000) 5159. 



21 



G. P. Beretta, J. Math. Phys. 25, (1984) 1507. 

K. Husimi, Proc. Phys. Math. Soc. Japan 22, (1940) 264. 

H. Fu and A. I. Solomom, J. Mod. Opt. 49, (2002) 259. 

Thomas Mainiero and Mason A. Porter, |arXiv:nlin/0702025fr 3 [nlin.CD] 20 Sep (2007). 

J. Rogel-Salazar, S. Choi, G.H.C. New, and K. Burnett, Phys. Rev. A 65, (2002) 023601. 

W. P. Reinhardt and C. W. Clark, J. Phys. B. 30, (1997) L785. 

L. M. Duan, A. Srensen, J. I. Cirac, P. Zoller, Phys. Rev. Lett. 85, (2000) 3991. 

A. Srensen, L.M. Duan , J. I. Cirac, P. Zoller, Nature 409, (2001) 63. 

M. S. Ateto, , arXiv: 0911.4240v2 [quant-ph], accepted in Int. J. Quant. Inf.. 

Amy Peng, Mattias Johnsson, W. P. Bowen, P. K. Lam, H. -A. Bachor and J. J. Hope, Phys. Rev. A 71, (2005) 033809. 

Jarorm'r Fiurasek, and Nicolas J. Cerf, Phys. Rev. Lett. 93, (2004) 063601. 

J. A. Vaccaro and A. Orlowski, Phys. Rev. A 51, (1995) 4172. 

A. Orlowski, H. Paul and G. Kastelewicz, Phys. Rev. A 52, (1995) 1621. 

I. Jex and A. Orlowski, J. Mod. Opt. 41, (1994) 2301. 

Andrew P. Hines, Ross H. McKenzie, and G.J. Milburn, |arXiv:quant-ph/0308165| ^4 28 Feb (2005) 



